Intrinsic rotation driven by non-Maxwellian equilibria in tokamak plasmas 
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The effect of small deviations from a Maxwellian equilibrium on turbulent momentum transport 
in tokamak plasmas is considered. These non-Maxwellian features, arising from diamagnetic ef- 
fects, introduce a strong dependence of the radial flux of co-current toroidal angular momentum on 
coUisionality: As the plasma goes from nearly coUisionless to weakly coUisional, the flux reverses 
direction from radially inward to outward. This indicates a collisionality-dependent transition from 
peaked to hollow rotation profiles, consistent with experimental observations of intrinsic rotation. 



Introduction. Observational evidence from magnetic 
confinement fusion experiments indicates that axisym- 
metric toroidal plasmas (tokamaks) that are initially sta- 
tionary develop differential toroidal rotation even in the 
absence of external momentum sources [IH5]. This 'in- 
trinsic' rotation can depend sensitively on plasma density 
and current, with relatively small variations reversing the 
rotation direction from co- to counter-current [21 |BHS]. 
Conservation of angular momentum dictates that the in- 
trinsic rotation is determined by momentum transport 
within the plasma. Since turbulence is the dominant 
transport mechanism in fusion plasmas [lOj , one must un- 
derstand turbulent momentum transport to understand 
intrinsic rotation. 

For the up-down symmetric magnetic equilibria used 
in most experiments, the turbulent momentum transport 
for a non-rotating plasma can be shown to be identically 
zero pTHIS] unless one retains formally small effects that 
are usually neglected in analysis. A self-consistent, first- 
principles theory has been formulated that includes these 
effects [131 [H]. Of these effects, only radial variation of 
plasma profile gradients [TBHTO] and slow variation of tur- 
bulence fluctuations along the mean magnetic field [20] 
have been studied, and these studies have not led to a 
theory that explains the key dependences of intrinsic ro- 
tation in the core of tokamaks. 

In this Letter we consider the novel effect of small 
deviations from an equilibrium Maxwellian distribution 
of particle velocities on turbulent momentum transport. 
These deviations arise naturally due to diamagnetic ef- 
fects in plasmas with curved magnetic fields and pressure 
gradients EI] . They vary strongly with quantities such as 
coUisionality, plasma current, and the equilibrium density 
and temperature gradients in the plasma. We show using 
direct numerical simulations that these non-Maxwellian 
features, though small, introduce significant new depen- 
dences to the turbulent momentum transport. We dis- 
cuss the physical origins of the dependences and possible 
implications for tokamak experiments. 

Momentum transport model. Tokamak plasma dy- 
namics typically consist of low amplitude, small scale 



turbulent fluctuations on top of a slowly evolving macro- 
scopic equilibrium. It is thus natural to employ a mean 
fleld theory in which the particle distribution function, 
/, is decomposed into equilibrium, F, and fluctuating, 
Sf, components. The fluctuations are low frequency, w, 
relative to the ion Larmor frequency, $7, and anisotropic 
with respect to the equilibrium magnetic field, with char- 
acteristic scales of the system size, L, along the field and 
the ion Larmor radius, p, across the field. Expanding 
/ = /o + /i + ..., employing the smallness parameter 
p* = p/L ~ uj/Q ~ Sf/F ~ fj+i/fj <^ 1, and averag- 
ing over the fast Larmor motion and over the fluctuation 
space-time scales, one obtains a coupled set of multiscale 
gyrokinetic equations for the fluctuation and equilibrium 
dynamics [221455] . 

Typically only the lowest order system of equations for 
Sf is considered. However, these equations have been 
shown to possess a symmetry that prohibits momen- 
tum transport in a non-rotating plasma [TTlllSj . Con- 
sequently, we include in our analysis higher order effects 
arising from corrections to the lowest order (Maxwellian) 
equilibrium [TH [IS]. We limit our analysis to these 
non-Maxwellian corrections because they are known to 
depend sensitively on plasma coUisionality and current, 
which are key parameters controlling intrinsic rotation 
in experiments [31 |5HS] • We further simplify analysis by 
considering only electrostatic fluctuations and by per- 
forming the subsidiary expansions p* ^ z^* ^ 1 and 
p* ^ Bg/B <^ 1, where B is the magnitude of the 
equilibrium magnetic field, Bg is the magnitude of the 
poloidal component, v^, == uaqR/vu, vu is the ion-ion 
collision frequency, vu = ^j2Ti/m.i is the ion thermal 
speed, Ti is the equilibrium ion temperature, mi is the 
ion mass, g is a measure of the pitch of the magnetic field 
lines called the safety factor, and R is the major radius 
of the torus. These are good expansion parameters in 
typical fusion plasmas. 

Using (R, e, p) variables, with R the position of the 
center of a particle's Larmor motion, e — mv^ /2 the par- 
ticle's kinetic energy, p — mv\/2B the particle's mag- 
netic moment, v the particle's speed, and _L indicating 



the component perpendicular to the magnetic field, the 
resulting equation for the fluctuation dynamics is 



( - rUsRvH 



(1) 



+ Z,ev||.Vl>^+a, 

where g — {5fi + (5/2), v is particle velocity, the sub- 
scripts II and _L denote the components along and across 
the equilibrium magnetic field, Ze is particle charge, 
(p — 54)1 + 5(t)2 and $ = $0 + ^1 are fluctuating 
and equilibrium electrostatic potentials, F — Fq + Fi, 
D /Dt ~ d/dt + v^; • Vj^, (.) is an average over Larmor 
radius at fixed R, v^j^ is the drift velocity due to the 
Coriolis effect and due to curvature and inhomogeneity 
in the equilibrium magnetic field, (5v£; — {c/B)h x W±(p 
and ve = {c/B)h x Vl> are E x B drift velocities, 
b is the unit vector directed along the magnetic field, 
uji^^E — ~{c/RBg){d^o/dr) is the toroidal rotation fre- 
quency due to E X B flow, c is the speed of light, the 
subscript s denotes species, and Cg describes the effect 
of Coulomb collisions on species s. 

Tokamak plasmas are sufficiently coUisional that the 
distribution of particle velocities is close to Maxwellian; 
i.e., fo~Ff)^ Fm, with Fm a Maxwellian. Equilibrium 
deviations from Fm are determined by the drift kinetic 
equation [21,, 



V|| • WH^s + VA/s • VF, 



Ms 



Cs[H- 



IsJ 



(2) 



where ^1=^1 + Zc^iFm/T. Finally, the electrostatic 
potentials are obtained and the system closed by enforc- 
ing quasineutrality: 

^Z,y'd3v(^5s + ^((^)-^)FM.)=0, (3) 
Y^ZsJ d'v (hu - j^-^iFms) = 0. (4) 



With gs and ip determined by Eqs. (Illl-Q, the tur- 
bulent radial fluxes of energy, Q, and toroidal angular 
momentum, 11, are given by 

Qs = {sJfsSvE ■ Vr)^ , (5) 

n = ^ {m^RHfs (v • VC) 5xE ■ Vr)^ , (6) 



where Sf = g + Ze{{ip) — (p)Fm/T, ( is toroidal angle, 
and (a)^ = J dt J d'^r J d^v a / J dt J d^r is an integral 
over all velocity space and over a volume of width w 
{p <S^ w <^ L) and time interval At {Ra/vn ^ Af <C 
p^'^Ro/vti) encompassing several turbulence correlation 
lengths and times. 



TABLE I 


: Collisionality dependence 


3f iJ^Cd 




i^t 


0.003 


0.030 


0.059 


0.089 


0.148 


0.208 


0.297 


Ro^cd 

Vti 


0.091 


0.114 


0.127 


0.137 


0.153 


0.165 


0.180 


Rl diucd 
Vti dr 


-0.447 


-0.577 


-0.651 


-0.701 


-0.776 


-0.829 


-0.891 



Results and analysis. We obtain the correction, Fi, 
to the equilibrium Maxwellian and the corresponding 
electrostatic potential, $1, by solving Eqs. ^ and Q 
using the drift kinetic code NEO [2^. These quanti- 
ties are then input to the Sf gyrokinetic code GS2 |28j . 
which we have modifled to solve Eqs. (IT]) and p) in 
the presence of Fi and <i>i. To calculate the 'intrinsic' 
momentum flux that is present even for a non-rotating 
plasma, we set the total toroidal angular momentum 
in a flux surface, which consists of diamagnetic and 
E X B contributions, to zero: J^s {{"^sR^^ ■ VC)/s) . = 
J2srf^sns{R'^)j^{oJc,E+^C.d) = 0, with Wf_d = 
J2s{'^sR'^i-y-'^C)Fis)^/J2s^sns{R'^)jy the diamag- 
netic contribution to the toroidal rotation frequency and 
n the number density. The non-zero E x B rotation 
needed to cancel the diamagnetic rotation breaks the 
symmetry of the lowest order gyrokinetic equation and 
thus contributes to momentum transport, as do the non- 
Maxwellian equilibrium corrections we have included. 

We consider a simple magnetic equilibrium with con- 
centric circular flux surfaces known as the Cyclone Base 
Case [2§\ . which has been benchmarked extensively in 
the fusion community. The equilibrium is fully specified 
by the Miller model 30], with q = lA, s = dlnq/dlnr = 
0.8, e = r/Ro = 0.18, Rq/L^ = 2.2, and Rq/Lt = 6.9, 
where: r is the minor radius at the constant pressure 
surface, or flux surface, of interest; Rq is the major ra- 
dius evaluated at r = 0; and L„ and Lt are the density 
and temperature gradient scale lengths for both ions and 
electrons. In order to obtain the gradient of Fi appearing 
in Eq. M , we must additionally specify the radial depen- 
dence of these quantities. For our base case, we choose 
Ro/Lm Rq/Lt, and dq/dr to be constant in radius. 

With these base case parameters specifled, we conduct 
a series of simulations with kinetic electrons and deu- 
terium ions, varying l'.^, and k == Rgd^ hiT/dr^ about the 
baseline value of i'^ = 0.003, k = 0. Our GS2 simula- 
tions use 32 grid points in the coordinate parallel to the 
magnetic field (the poloidal angle), 12 grid points in e, 
37 grid points in A = fi/e, and 128 and 22 Fourier modes 
in the radial and binormal coordinates, respectively. The 
box size in both the radial and binormal coordinates is 
approximately 125^.^. 

The resulting H/Qi values as a function of 1^^ are 
shown in Fig. IT] We normalize H by Qi, which is al- 
ways positive, to remove any dependence of overall tur- 
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FIG. 1: Ratio of radial fiuxes of ion toroidal angular mo- 
mentum, n, and energy, Qi, vs. normalized ion-ion collision 
frequency, Vt- 



bulence amplitude on collisionality. Note that Fi, and 
thus UJ^.E — —^Q,d^ varies with coUisionality, as indicated 
in Table llj For nearly coUisionless plasmas, II/Qi is neg- 
ative, indicating a radially inward flux of co-current an- 
gular momentum that would contribute to a centrally 
peaked rotation profile. The ratio II/Qj increases with 
I'*, passing through zero and becoming positive when 
I/* ~ e'^/^. For v^ > e"!"^ , the radially outward flux of co- 
current angular momentum would contribute to a hollow 
rotation profile. 

In Fig.[2J we show results from a series of simulations in 
which we independently set the ExB rotation (including 
its derivative) and the diamagnetic effects, represented 
by Fi, to zero. These are given by the blue and red 
curves, respectively. We see that the ExB rotation 
causes an inward momentum flux, with Tl/Qi increasing 
in magnitude with v^. The non-Maxwellian correction 
Fi gives a H/Qi that goes from slightly negative to large 
and positive as v^, is increased. A partial cancellation 
between these effects gives the actual H/Qi. 

To explore in more detail the origin of the sign reversal 
of n/Qi, it is convenient to express the ion energy flux 
in the diffusive form Qi ~ —XidTi/dr, and to decompose 
the momentum flux into diffusive, advective, and other 
pieces: 
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dr '^"■'^ dr '""■^J (7) 

- mRo {oJt:,dPd + uji^^ePe) + Ilothor, 



where x<p,d and Pd are diffusion and advection (commonly 
called 'pinch') coefficients, respectively, for the diamag- 
netic rotation, and X4>,e and Pe play the same roles for 
the ExB rotation. The quantity Ilothor accounts for 
all other sources of 11 that arise due to Fi(cL'^^(r) — 
i-^CEif) = 0); e.g., the equilibrium parallel heat flow and 
other higher order velocity moments of Fi will contribute 
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FIG. 2: Ratio of radial fluxes of ion toroidal angular mo- 
mentum, n, and energy, Qi, vs. normalized ion-ion collision 
frequency, i/, for: the base simulations (black), simulations 
with no E X B rotation to balance the diamagnetic rotation 
(blue), simulations with no correction, Fi, to the Maxwellian 
equilbrium (red), and simulations with Ilother = (green). 



to Ilothcr- Using the fact that ujq^e 
rotating plasma, we have 
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where X0,cff = X4>.d - X<t>.E and Pcff = Pd- Pe- 

Changing i>^ can alter H/Qi in multiple ways. First, 
the turbulent advection and diffusion coefficients, Pcff/Xi 
and X4>,cs/Xi7 can be modified either directly by colli- 
sions or indirectly by the :/* -dependent rotation and ro- 
tation gradient. By independently varying i/*, w^^^;, and 
dujQE/clr in GS2 turbulence simulations with fixed Fi 
and $1, we found that such modifications of the ad- 
vection and diffusion coefficients were minor. Further- 
more, the turbulence type, characterized by the dom- 
inant linear instability mechanism, remained the same 
(ion-temperature-gradient driven) for all simulations. 

With x<i>,es/Xi and Pcs/Xi approximately independent 
of i^* , we see from Eq. (Is]) that the j/* dependence of (11 — 
nother)/Qi comes entirely from the change of ujc^^d and 
dujQ^d/dr with v^., given in Table ^ In order to calculate 
(n — Ilothcr )/Qj I we ran a series of simulations in which 
we used a modified Fi that was constrained to produce 
pure rotation so that Ilother = 0. The results are shown 
as the green curve in Fig. [2] We see that (11 — Ilother)/*?! 
is always positive and increases approximately linearly 
with dujQ,d/dr, as diffusion was found to dominate over 
advection in these cases. This indicates that equal and 
opposite diamagnetic and ExB rotations do not lead 
to a complete cancellation of momentum transport |31j . 
The increase in (11 — Ilother)/'?* with v.^, accounts for just 
over half of the total increase in li/Qi over the range of 
i^* we have considered. The rest of the increase, as well 
as the negative offset needed to give the sign reversal in 
li/Qi must come from IIothor/Qi- 



TABLE II: Temperature profile dependence of o;^., 
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-1.116 


-0.447 


-0.105 


0.223 


0.835 



To see how these resuhs may be modified for dif- 
ferent plasma profiles, we also conducted a series of 
simulations in which we fixed v^, — 0.003 and varied 
K = R^d"^ h\T / dr"^ . Since the calculation of Fi in NED 
depends on Rq/Lt, varying k affects dFi/dr but not Fi 
itself. Consequently, duj^,d/dr varies with n (see Table [Tl|) 
while uj(^^d itself remains fixed. The change in H/Qi with 
K is shown in Fig. [3J As was the case in the z^* study, 
the E X B and Fi contributions to H/Qi partially can- 
cel, though in this case each contribution independently 
changes sign with k. The net result is a relatively weak 
variation of Tl/Qi with no sign reversal. 

Discussion. The sign reversal of H/Qi shown in Fig.fl] 
suggests a transition from peaked to hollow rotation pro- 
files when z/* ^ e"^' ^. This is consistent with experimental 
results, which show such transitions at similar i^^ values 
when density (proportional to v^) is increased or current 
(inversely proportional to i^* ) is decreased [31 O 111 |32] . 
Furthermore, our observation that the normalized tur- 
bulence diffusion and advection coefficients vary only 
minimally during the transition agree with recent ex- 
perimental observations showing that the fundamental 
turbulence characteristics are unaltered as the rotation 
reverses direction 0. 

From Fig. [2] and the analysis following Eq. (|8| , it is 
evident that a combination of effects leads to the sign re- 
versal of H/Qi- However, the sign reversal fundamentally 
originates from the j/* dependence of Fi , which has been 
extensively studied and is the main concern of 'neoclas- 
sical' theory (see, e.g., [2TJ[33]). For v^, <^ e^^"^, known 
as the 'banana' regime, all particle orbits are collision- 
less. However, for e"^/^ <^ v^, <S^ I, known as the 'plateau' 
regime, low energy particles that are trapped in the equi- 
librium magnetic well become coUisional. For a plasma 
perfectly in the banana or plateau regimes, one can show 
that -Fi, and thus our 11/ Qi, becomes independent of 
f^ [m |33] . It is only when transitioning between these 
regimes that 11/ Qi varies with j/*. So, while different 
profiles of quantities such as density, temperature, and 
current may alter or eliminate the transitions with v^, 
discussed above, they can only occur for v^ ^ e"!"^ . 

Finally, we reiterate that in our analysis we retained 
small terms (namely the diamagnetic effects that give 
rise to departures from a Maxwellian equilibrium distri- 
bution) in the multiscale gyrokinetic expansion, while we 
neglected other terms (radial profile variation, certain ef- 
fects arising from the slow variation of fiuctuations along 
the magnetic field, etc.) that may be of the same size. 
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FIG. 3: Ratio of radial fiuxes of ion toroidal angular mo- 
mentum, n,;, and energy, Qi, vs. normalized second deriva- 
tive of the logarithmic ion temperature for: the base simu- 
lations (black) and for simulations with no E x B flow to 
balance the diamagnetic flow (red) and no correction, Fi, to 
the Maxwellian equilibrium (blue). 



There are two justifications for this. First, if the fluc- 
tuation amplitudes and scales do not vary strongly with 
Bg/B, then the diamagnetic effects considered here dom- 
inate so that our model is fully self-consistent [TH [TS] . 
Second, the small effects we have neglected are not ex- 
pected to have a particularly strong dependence on coUi- 
sionality. Thus, while inclusion of these effects may pro- 
vide an offset to the momentum transport, we do not 
expect them to modify the variation of Tl/Qi with v^ 
presented here. 
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